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ABSTRACT 

We present a three-dimensional (3-D) dynamical model for the broad [Fe ill] emission ob- 
served in 1] Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectro- 
graph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) 
simulations of rj Car's binary colliding winds. Radiative transfer codes are used to generate 
synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases 
and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclina- 
tion i, the PA 6 that the orbital plane projection of the line-of-sight makes with the apastron 
side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first 
time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously re- 
produce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38°, 
and the temporal variations in emission seen at negative slit PAs, the binary needs to have an 
i « 130° to 145°, 9 Ri -15° to +30°, and an orbital axis projected on the sky at a PA « 302° 
to 327° east of north. This represents a system with an orbital axis that is closely aligned with 
the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, thus orbits 
clockwise on the sky and is on the observer's side of the system at apastron. This orientation 
has important implications for theories for the formation of the Homunculus and helps lay the 
groundwork for orbital modeling to determine the stellar masses. 
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' 1 INTRODUCTION 

?7 Carinae is the most luminous, evolved stellar object that can be 
closely studied (Davidson & Humphreys 1997). Its immense lu- 
minosity (L > 5 X 10'' 1/0, ICoxet al. 1995) and relative prox- 
imity (D = 2.3 ±0.1 kpc. ISmith 2006b) make it possible to 
test and constrain specific theoretical mo dels of extremely massive 
stars (> 6 Mq) using high-quality data jHillier et al.ll200ll , l200d : 

ISmith & O wocki 2006). rj Car is thus one of the most intensely 
observed stellar systems in the Galaxy, having been the focus of 
numerous ground- and space- based obse rving campaign s at mul- 
tiple wavelengths fro m radio jPuncan & W hite 2003; W hite et al.1 

l2005t) to gamma-rays jXavani et al.ll2009l : lAbdo et alj|201ol) ^ 



* Based on observations made with the NASA/ESA Hubble Space Tele- 
scope. Support for programs 7302, 8036, 8483, 8619, 9083, 9337, 9420, 
9973, 10957 and 11273 was provided by NASA directly to the STIS Sci- 
ence Team and through grants from the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc., under NASA contract NAS 5-26555. 
f Email: tmadura@mpifr-bonn.mpg.de 
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Unfortunately, the dusty Homunculus nebula that formed dur- 
ing -q Car's "Great Eruption" in the 1840s enshrouds the s ystem, 
complicating direct observations of the central stellar source ismithl 
1 200^) . Nevertheless, ground- and space-based, multi-wavelength 
observations obtained over the past two decades strongly indicate 
that rj Car is a highly eccentric (e ~ 0.9) colliding wind binary 
(CWB) with a 5.54-year orbital period (.Damineli 1996: Feast et al. 



| 200l| : JPunc an & Wh it3l2003l: IWhitelock et alj|2004 ISmith et al 



2004"; "Corcoran 2005: Vemer et al.'2005':'yan Ge nderen et al.l2006l. 
Damineli et al. 2008a b; Fernandez-Laiijs et alj I2OIOI : ICorcoranI 
I2OIIL hereafter CI 1). The consensus view is th at the primary star , 
r?A, is a Luminous Blue Vari a ble (LBV) tHillier & AllenI 1 19921 : 
iDavidson & Humphreys! 1 1991 iHillier et al.1 1200 iL l2006h Radia- 
tive transfer modeling of Hubble Space Telescope/Space Tele- 
scope Imaging Spectrograph (HST/STIS) spatially-resolved spec- 
troscopic observations suggests that riA has a current mass > 
90 Mq, and a current-day stellar wind with a mass-loss rate of 
y lO'^Mp yr"^ and t erminal speed of ~ 500 - 600 kms"^ 
jHillier et alj200 l'.'2006'. hereafter HOI, H06). The companion star, 
?7b, continues to evade direct detection since tja dwarfs its emission 
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at most wavelengths (Damineli et al. 2000; HOI; H06; Smith et al. 
2004; Cll). As a result, ?7b's stellar parameters, evolutionary state, 
orbit, and influence on the evolution of 77A are poorly known. 

Currently, the best constraints on the stellar properties of 
r(B come from photoionization modeling of the "Weigelt blobs", 
dense, slow-moving ejecta in the vicinity of the binary system 
dWeigelt & Ebersbergej[T983) . I Vemer et all ( I2OO5I) initially char- 
act erized r?B as a mid-O to WN supergiant. More recent work 
bv lMehneretZ I Jioia^ tightly constrains the effective tempera- 
ture of ?7B (Teff ~ 37, 000 - 43, 000 K), but not its luminosity 
(log Z//L0 ~ 5 — 6), resulting in a larger range of allowed stellar 
parameters. 

Constraints on the wind parameters of rye come from ex- 
tended X-ray monitoring by the Rossi X-ray Timing Explorer 
(RXT E ), Chandra, XMM, and Suzaku sa t ellites alshibashi et al. 

lb 



19991; ICorcoranI |2005| ; iHamaguchi etd] |2007| ; fcorcoran et al] 



201(1) . The periodic nature, minimum around periastron, and hard- 



ness (up to 10 keV) of 77 Car's RXTE light curve are all character- 
istics of a highly eccentric CWB, the variable X-ray emission aris- 
ing in a wind- wind collision (WWC) zone formed between the two 
stars ('Pittard & Corcoran 2002; Okazaki et aLll2008l ; IParkin et al] 
[2009, 2011; Corcoran et al. 2010, hereafter PC02, O08, P09, Pll, 
and CIO, respectively). The hardness of the X-rays requires that 773 
have a high wind terminal speed of ~ 3000 km s"'^, and detailed 
modeling of the momentum balance between the two shock fronts 
suggests a mass-loss rate of ~ lO^^M© yr"^ (PC02; O08; P09; 
Pll). 

Proper numerical modeling of rj Car's WWC remains a chal- 
lenge, mainly because it requires a full three-dimensional (3-D) 
treatment since orbital motion, especially during periastron, can be 
important or even dominant, affecting the shape and dynamics of 
the WWC region (O08; P09; Pll). Most hydrodynamical simula- 
tions of rj Car have been two-dimensi onal (2-D), neglect ing the 
effects of orbital motion for simplicity iPittardll998l.l200(i PC02; 
Henley 2005). Until very recently (O08; Pll), fully 3-D simula- 
tions were computationally impractical. 

When it comes to reproducing observational diagnostics of 
r] Car from hydrodynamical simulations of its colliding winds, the 
focus has been almost exclusively on X-rays. One drawback of this 
is that models for 77 Car at other wavelengths have been mostly phe- 
nomenological, offering only qualitative explanations for the vast 
array of complicated observations. This makes it difficult to test and 
refine models, or quantitatively constrain the physical parameters of 
the system. An excellent example is the lack of consensus regard- 
ing rj Car's orbital orientation. The majority favor an orbit in which 
77B is behind j^a during periastron (Damineli 1996; PC02; Corco- 
ran 2005; Hamaguchi et al. 2007; Nielsen et al. 2007a; Damineli et 
al. 2008b; Henley et al. 2008; O08; P09; Pll; Groh et al. 2010b; 
Richardson e t al. 2010). But others place t^b on the near side of r]A 
at periastron (jpalceta-Goncalves et al.l2005l;lAbraham et al 2005t ; 



Abraham & Falceta-Goncalves|l2007l ; iKashi & Sokerl 120071 12008 ; 
Falceta-Goncalves & AbrahamI 20090. Some claim that the orienta- 



tion is not well established at all jMehner et al. 201 1). 

A more fundamental drawback of focusing purely on X-rays 
or other spatially-unresolved data is that any derived orbit is am- 
biguous with respect to its absolute orientation on the sky; any or- 
bital orientation derived solely from fitting the X-ray data can be 
rotated on the sky about the observer's line-of-sight and still match 
the observations. A degeneracy also exists in the orbital inclination, 
with models that assume i ~ 45° and i ~ 135° both capable of fit- 
ting the observed RXTE light curve. However, high-resolution spa- 
tial information is not enough; moderate-resolution spectral data is 



needed to determine the velocity structure of the emitting gas, and 
thus fully constrain the orbital orientation. As a result, using X-ray 
data alone, it is impossible to determine the 3-D alignment or mis- 
alignment of the orbital axis with the Homunculus polar axis, or 
the direction of the orbit. It is commonly assumed that the orbital 
axis is aligned with the Homunculus polar axis, but to date, neither 
modeling nor observations have unambiguously demonstrated that 
this is the case. Knowing the orbit's orientation is key to constrain- 
ing theories for the cause of the Great Eruption and formation of the 
Homunculu s (i.e. single star versus binary interaction/merger sce- 



narios.'lben"l999';'Dwarkadas & Owocki"2002';' Smith et alj200"33 ; 
Soker 2007; Smith 2009, 201 1 ; Smith & Frew 20Tl|). Furthermore, 
a precise set of orbital parameters would lay the groundwork for 
orbital modeling to determine the stellar masses. 

In a recent attempt to characterize rj Car's interacting 
winds, iGull et alj ( I2OO9I , hereafter G09) presented an analysis of 
HST /STIS observations taken between 1998 and 2004, identifying 
spatially-extended (up to 0.8"), velocity-resolved forbidden emis- 
sion lines from low- and high-ionizatiorlj species. HST/STIS imag- 
ing spectroscopy is ideal for moving beyond X-ray signatures as it 
uniquely provides the spatial (~ 0.1") and spectral (R ~ 8, 000) 
resolution necessary to separate the spectra of 77 Car's central 
source and nearby circumstellar ejecta from those of the Homuncu- 
lus and other surrounding material. Moreover, this high-ionization 
forbidden emission is phase-locked, which strongly suggests it is 
regulated by the orbital motion of tjb ■ Observed spectro-images are 
also highly dependent on the P/^ of the HST/STIS slit, implying 
that they contain valuable geometrical information. This opens up 
the possibility of obtaining 77B 's orbital, stellar, and wind parame- 
ters through proper modeling of the extended forbidden line emis- 
sion. 

This paper presents a detailed 3-D dynamical model for the 
broad components of the high-ionization forbidden line emission 
observed in 77 Car using the HST/STIS. Our model is based on the 
results of full 3-D Smoothed Particle Hydrodynamics (SPH) simu- 
lations of 77 Car's colliding winds (§ 3.2 and 4.1). Radiative trans- 
fer ca lculations p erformed with a modified version of the SPLASH 
code ( |Pricell2007l) are used with IDL routines to generate synthetic 
spectro-images of [Fe III] emission line structures at various orbital 
phases and STIS slit PAs (§ 3.4). Through a parameter study that 
varies the orbital inclination, the angle that the line-of-sight makes 
with the apastron side of the semi-major axis, and the PA on the 
sky of the projected orbital axis, we are able to tightly constrain, 
for the first time, the absolute (3-D) orientation and direction of 
77 Car's orbit (§ 4), showing in particular that the orbital axis is 
closely aligned in 3-D with the inferred polar axis of the Homuncu- 
lus nebula. A discussion of the results and this derived orientation 
are in § 5. § 6 summarizes our conclusions and outlines the direc- 
tion of future work. We begin (§ 2) with a brief summary of the 
observations used in this paper. 



2 THE OBSERVATIONS 

This paper uses the same HST/STIS data described in G09 (see their 
Table 1). Extracted portions of observations recorded from March 
1998 to March 2004 with the STIS CCD (0.0507" pixeP^ scale) 
with medium dispersion gratings {R ~ 8, 000), in combination 



^ Low- and high-ionization refer here to atomic species with ionization 
potentials (IPs) below and above the IP of hydrogen, 13.6 eV. 
^ Position Angle, measured in degrees from north to east. 
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Figure 1. Examples of the HST/STIS spatially-resolved [Fe III] spectra. Top Left: An //ST/ Advanced Camera for Surveys (ACS) High Resolution Camera 
(HRC) image of the Homunculus and r] Car with a 20" tield of view. Bottom Left: Enlarged central 2" field of view with the 52" X 0.1" slit positioned 
at PA = +69° and directional labels (NW = northwest, etc.) included for reference in the text. Spectro-images in the middle and right-hand columns were 
recorded from the central 2" X 0.1" portion of the slit, as drawn between labels I and II. Top Center: Spatially-resolved line profile of [Felll] A4659 in 
original form, recorded at PA = +69° in 2002 July (</> = 0.82). Top Right: The same spectro-image with the continuum subtracted on a spatial row-by-row 
basis. The velocity scale, referenced to the vacuum rest wavelength of the line, is ±600 km s^^, and the color level is proportional to the square root of the 
intensity. Bottom Center: Spectro-image of [Fell] A4815. shifted by 400 km s~^, used as a template to determine the spatial and spectral ranges where the 
[Fe III] image is contaminated by [Fe II] A4666 emission. Bottom Right: Same specfi'o-image as the top right panel, but with masks that remove the narrow 
line contamination and the contamination to the red from nearby [Fe II] A4666 (see text). 



with the 52" x 0.1" slit from 1,640 to 10,100 A, were used to sam- 
ple the spectrum of the central core of Car. Phases of the obser- 
vations^, arerelative to the RXTE X-ray minimum at 1997.9604 
( ICorcoraij|2005h : JDobs = JD 2450799.792 + 2024.0 x <;/>. Solar 
panel orientation requirements constrained the STIS slit PAs. The 
reduced STIS CCD spectra available through the STScI archives 
( |http : / / archive .stsci.edu /prepds/etacar/) were 
used. Compass directions (NW = north-west, etc.) describe the spa- 
tial extent of the emission. All wavelengths are in vacuum and ve- 
locities are heliocentric. For details, see G09 and references therein. 

Figure [T] shows examples of resolved, broad emission from 
[Fe III] A4659. In the top center panel, faint spatially- and velocity- 
resolved emission can be seen against bright nebular emission and 
dust-scattered stellar radiation. Strong, narrow line emission from 
the SW of rj Car is also seen. The bright continuum close to the 
wavelength of the wind line of interest is subtracted on a spatial 
row-by-row basis in order to isolate the fainter broad forbidden 
line emission (top right panel of Figure [TJ. All spectro-images in 
this paper have been processed identically using a portion of the 
spectrum with no bright narrow- or broad-line contamination. Row- 
by-row subtraction across the stellar position is less successful due 
to insufficient correction in the data reduction for small tilts of the 
spectrum on the CCD (G09). 

The bottom right panel of Figure[T]contains the same spectro- 
image of [Fe III] A4659 as the top right panel, but with masks that 
remove narrow line contamination from the Weigelt blobs, and con- 



tamination to the red due to the nearby [Fell] A4666 wind line. 
The masks are intended to remove "distracting" emission features 
not relevant to the work in this paper, displaying only the observed 
broad structures that are the focus of the modeling (outlined in 
white). To determine the spatial and spectral ranges where [Fell] 
A4666 contaminates the [Felll] A4659 image, we use as a tem- 
plate the uncontaminated spectro-image of the bright [Fell] A4815 
line (bottom center panel of Figure[Tll. We also compare the [Fe III] 
A4659 images to those of the [Fe III] A4702 and [N II] A5756 lines, 
which form in nearly identical conditions as they have very similar 
critical densities and IPs (Table 3 of G09). A detailed discussion 
of the observations and masking procedure for all spectro-images 
modeled in this paper is in Appendix A (available in the online 
version of this paper - see Supporting Information). 



2.1 The Key Observational Constraints JVIodeled 

Three key features observed in spectro-images of the high- 
ionization forbidden line emission can be used to constrain the 3-D 
orientation of rj Car's orbit. Each provides important clues about 
the nature and orbital variation of rj Car's interacting winds. We 
briefly summarize these below. 
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2.1.1 Constraint 1: Emission Arcs at Slit PA = +38°, </> = 0.976 

Figure |2] displays observed spectro-images of [Felll] A4659, 
recorded at slit PA = +38° on 2003 May 5 (0 = 0.976), sliow- 
ing spatially-extended (up to ~ 0.35") emission in the form of 
very distinct, nearly complete arcs that are entirely blue-shifted, 
up to ~ —475 km s^^. Several weaker [Felll] lines nearby in 
the spectrum, while blended, show no evidence of a spatially- 
extended (> 0.1") red component, and no other high-ionization 
forbidden lines show an extended red component at this position 
(G09). To illustrate this, the bottom right panel of Figure[2]contains 
a continuum-subtracted spectro-image of [N II] A5756, taken at the 
same orbital phase and slit PA. The IP (14.5 eV) and critical density 
(3 X 10^ cm"^) of [N II] A5756 are very similar to that of [Fe III] 
A4659 (16.2 eV and ~ 10^ cm^'^), meaning that their spatially- 
extended emission forms in nearly identical regions and via the 
same physical mechanism. The [N II] A5756 line is strong and lo- 
cated in a spectral region with very few other lines. Unlike [Fe III] 
A4659, there is no contamination to the red of [N II] A5756, mak- 
ing it ideal for demonstrating that there is no spatially-extended, 
red-shifted, broad high-ionization forbidden line emission. 

Of particular importance are the asymmetric shape and inten- 
sity of these emission arcs. Both arcs extend to nearly the same 
spatial distance from the central core, but the upper arc is notice- 
ably dimmer and stretches farther to the blue, to ~ —475 km s^^. 
The lower arc is brighter, but only extends to ~ —400 km s"'^ . Ev- 
ery high-ionization forbidden line observed at this orbital phase and 
slit PA exhibits these asymmetries in intensity and velocity (figures 
6 — 9 of G09), indicating that they are independent of the IP and 
critical density of the line and intrinsic to the shape and distribution 
of the photoionized extended wind material in which the lines form. 
As such, the shape and asymmetry of the blue-shifted emission arcs 
may be used to help constrain the orbital orientation. 

2.1.2 Constraint 2: Variations with Orbital Phase at Constant 
Slit PA = -28° 

Let us focus next on shifts of the high-ionization forbidden emis- 
sion with orbital phase. Figure 11 of G09 shows sets of six spec- 
tra of [Fe III] A4659 and A4702 recorded at select phases at slit 
PA = -28°. During the two minima at = 0.045 and 1.040 (rows 
1 and 5 of figure 1 1 of G09), the broad [Fe III] emission is absent. 
By = 0.213 (row 2) and <^ = 1.122 (row 6, ~ 8 months after the 
X-ray minimum), it strongly reappears. Images at </!> — 0.407 are 
similar to those at (p = 0.213 and 1.122. By <j) = 0.952, emission 
above ~ 150 km disappears. This phase dependence in the 
[Felll] emission, especially the disappearance during periastron, 
strongly indicates that it is tied to the orbital motion of ijb . 

Another important observational feature is the presence of 
spatially-extended red-shifted emission to the NW, in addition to 
extended blue-shifted emission to the SE, at phases far from pe- 
riastron. This is in contrast to the entirely blue-shifted arcs seen 
for slit PA = +38°, indicating that the STIS slit at PA = -28° 
is sampling different emitting portions of rj Car's extended wind 
structures, even at similar orbital phases. 

2.1.3 Constraint 3: Doppler Shift Correlations with Orbital 
Phase and Slit PA 

The final set of observations focuses on changes with slit PA lead- 
ing up to periastron. Figures 12 and 13 of G09 illustrate the be- 
havior of the [Fe III] emission for various slit PAs at select phases. 




Velocity (km/s) Velocity (km/s) 

Figure 2. Spectro-images recorded 2003 May 5 (fli = 0.976). Top Left 
Panel: HST/ ACS HRC image of the central 2" of -q Car with the 2" long 
portion of the STIS slit at PA = +38° overlaid. The locations of Weigelt 
blobs B, C, and D are indicated, as is the direction of north. Top Right Panel: 
Continuum-subtracted spectro-image of [Fe III] A4659. Emission appears 
as a pair of completely blue-shifted arcs. Color is proportional to the squai'e 
root of the intensity and the velocity scale is ±600 km s~ ^ . Weak emission 
to the red is entirely due to contamination by [Fe 11] A4666, while weak, 
nan'ow emission is centered at ~ —40 km s^^. Bottom Left Panel: Same 
spectro-image, but with masks that remove the line contamination. Bottom 
Right Panel: Continuum-subtracted spectro-image of [N II] A5756 included 
to demonstrate that at this phase and STIS PA, the high-ionization forbidden 
lines show no evidence of broad, spatially-extended, red-shifted emission. 

The key point to take away from this particular data set is that away 
from periastron, the spatially-extended [Fe III] emission is almost 
entirely blue-shifted for positive slit PAs = +22° to +70°, but is 
partially red-shifted for negative PAs, most notably PA = —82°. 



3 A 3-D DYNAMICAL MODEL FOR THE BROAD 
HIGH-IONIZATION FORBIDDEN LINE EMISSION 

3.1 Why Model [Fe ill] A4659 Emission? 

Table 3 of G09 lists the various forbidden emission lines observed 
in rj Car. However, there are several important reasons to focus par- 
ticular attention on [Fe III] A4659 emission. 

First, the physical mechanism for the formation of forbidden 
lines is well understood (Appendix B, online version). Moreover, 
emission from forbidden lines is optically thin. Thus, forbidden line 
radiation can escape from a nebula much more easily than radiation 
from an optically thick resonance line, which is emitted and reab- 
sorbed many times (Hartman 2003). This is of great value when 
studying an object so enshrouded by circumstellar material since 
one does not have to model complicated radiative transfer effects. 

[Felll] is considered rather than a lower ionization line of 
[Fe II] because the bulk of the [Fe III] emission arises in regions that 
are directly photoionized by rjB (Verner et al. 2005; G09; Mehner et 
al. 2010). No intrinsic [Felll] emission is expected from rjA- This 
is based on detailed theoretical models by HOI and H06, which 
show that in tja's envelope, Fe'^"'" only exists in regions where the 
electron density is two-to-four orders of magnitude higher than the 
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critical density of the [Felll] A4659 line. Since Fe" needs 16.2 eV 
radiation (or collisions) to reach the Fe'^^ state, most of the [Fe III] 
emission arises in areas directly photoionized by rjB ■ 

In comparison, Fe" only needs 7.9 eV radiation or collisions to 
ionize to Fe^, which can form in the wind of rjA (without r^B S in- 
fluence) and/or near the wind-wind interaction regions in areas ex- 
cited by mid-UV radiation filtered by 77a's dense wind. Collisions 
and photoexcitation to upper Fe^ energy levels can also populate 
many metastable levels. As a result, [Fell] emission is far more 
complicated and originates from lower excitation, lower density re- 
gions on much larger spatial scales (HOI; H06; G09; Mehner et al. 
2010). Because of this complexity, the modeling of [Fe II] emission 
is deferred to future work. 

The IP of 16.2 eV required for Fe^+ also means that the 
[Fe III] emission forms in regions where hydrogen is ionized, but 
helium is still neutral. In contrast, both Ar'^+ and Ne^+ have IPs 
greater than 24.6 eV. Therefore, [Ar III] and [Ne III] emission arise 
in areas where helium is singly ionized. Modeling [Fe III] is thus 
more straightforward as one does not have to worry about the two 
different types of ionization structure possible (one due to H and 
one due to He), which depend on the spectrum of ionizing radia- 
tion and the abundance of helium ( Osterbrock 1989). 

[Fe III] A4659 is chosen over the similar emission line of [N II] 
A5756 in order to avoid potential complications due to intrinsic 
[Nil] emission in the extended wind of rjA- [Nil] emission can 
form in a broad zone between 100 and 1000/?, (~ 30 - 300 AU) 
of TjA and is sensitive to its wind temperature and mass-loss rate 
(HOI). By focusing on [Felll], one avoids having to include any 
intrinsic forbidden line emission from tja- 

The modeling in this paper focuses solely on the broad [Fe III] 
emission features that are thought to arise in the dense, mod- 
erate velocity (~ 100 — 600 kms^^) extended primary wind 
and WWC regions. The much narrower (< 50 kms"^) emis- 
sion features that form in the Weigelt blobs and other dense, 
slow-moving equatorial ejecta from rj Car's smaller eruption 
in the 1890s [W eigelt & Ebersberger 1986"; 'Davidson et al. 1995", 



Table 1. Stellar, Wind, and Orbital Parameters of the 3-D SPH Simulation 



I997I: llshibashi et ak .2003, ; Smith et ak 2004; Gull et al. 200ft ; 
Mehner et al.ll201ol) are not modeled. Any possible effects of local 



dust formation are neglected for simplicity as not enough informa- 
tion is available to realis tically include them at this time ( Williams 
l2008l ; lMehner et alj2OI0b . 

As all observed high-ionization forbidden lines show the same 
basic spatial and temporal features in their broad emission, with any 
differences in size or location attributable to differences in the IP 
and critical density of the specific line of interest, focusing solely 
on [Fe III] A4659 should not significantly bias the overall conclu- 
sions, which should extend to the other high-ionization forbidden 
lines as well. Values of the transition probability {A21 = 0.44 s~^), 
statistical weights (gi — g2 ~ 9), and line transition frequency 



(V21 



6.435 X 10^"^ s-i) for [Felll] A4659 all come from 



iNahar & PradhanI ( Il996h . Values for the collision strengths {fli2) 
come from Zhand ( 1996h . A sola r abundance of iron is also used, 
consistent with the works of HOI; IVemer et alj ( l2005h ; and H06. 



3.2 The 3-D SPH Code and General Problem Setup 

The numerical simulations in this paper were performed with the 
same 3-D SPH code used in O08. The stellar winds are modeled 
by an ensemble of gas particles that are continuously ejected with 
a given outward velocity at a radius just outside each star. For sim- 
plicity, both winds are taken to be adiabatic, with the same initial 
temperature (35, 000 K) at the stellar surfaces and coasting with- 



Parameter 


VA 


Vb 


Mass (Mq) 


90 


30 


Radius (Rq) 


90 


30 


Mass-Loss Rate {Mq yr~^) 


10" 


^ 10-5 


Wind Terminal Velocity (km s~^) 


500 


3000 


Orbital Period (days) 




2024 


Orbital Eccentricity e 




0.9 


Semi-major Axis Length a (AU) 




15.4 



out any net external forces, assuming that gravitational forces are 
effectively canceled by radiative driving terms (O08). The gas has 
negligible self-gravity and adiabatic cooling is included. 

In a standard xyz Cartesian coordinate system, the binary or- 
bit is set in the xy plane, with the origin at the system center-of- 
mass and the orbital major axis along the s-axis. The two stars 
orbit counter-clockwise when looking down on the orbital plane 
along the +z axis. Simulations are started with the stars at apas- 
tron and run for multiple consecutive orbits. By convention, t — 
(4> = 0) is defined to be at periastron passage. 

To match the ~ ±0.7" scale of the STIS observations a t the 
adopted distance of 2.3 kpc to the 77 Car system ( ISmithll2006bh . the 
outer simulation boundary is set at r = ±105a from the system 
center-of-mass, where a is the length of the orbital semi-major axis 
(a — 15.4 AU). Particles crossing this boundary are removed from 
the simulation. We note that the SPH formalism is ideally suited 
for such large-scale 3 -D simulations, as compared to grid-based 
hydrodynamics codes ( |Pricell2004l ; lMonaghan|[200^ . 

Table [T] summarizes the stellar, wind, and orbital parameters 
used in the modeling, which are consistent with those derived from 
the observations (PC02; HOI; H06; CI 1), with the exception of the 
wind temperature of rjA- The effect of the wind temperature on the 
dynamics of the high-velocity wind collision is negligible (O08). 
Note that the simulation adopts the higher mass- loss rate for rjA 
derived bv lDavidson et aI.I ( ll995h . lCox et alj|l995b . HOI, and H06, 
rather than the factor of four lower mass-loss rate assumed when 
analyzing X-ray signatures (PC02; O08; P09; PI 1; CI I). Effects of 
the adopted value of the primary mass-loss rate on the forbidden 
line emission will be the subject of a future paper^ 

The publicly available software SPLASH ( Pricel2007l) is used 
to visualize the 3-D SPH code output. SPLASH differs from other 
tools because it is designed to visualize SPH data using SPH algo- 
rithms. There are a number of benefits to using SPLASH, and the 
reader is referred to Price (2007) and the SPLASH userguide for 
these and discussions on the interpolation algorithms. 



3.3 Generation of Synthetic Slit Spectro-Images 

The 3-D SPH simulation provides the time-dependent, 3-D den- 
sity and temperature structure of 77 Car's interacting winds on spa- 
tial scales comparable to those of the HST/STIS observations. This 
forms the basis of our 3-D dynamical model. Unfortunately, it is 
currently not possible to perform full 3-D simulations of 77 Car in 
which the radiative transfer is properly coupled to the hydrodynam- 
ics ( Paardekooper 2010). Instead, the radiative transfer calculations 
here are performed as post-processing on the 3-D SPH simulation 
output. This should not strongly affect the results so long as the ma- 
terial that is photoionized by 773 responds nearly instantaneously to 
its UV flux. In other words, as long as the timescale for the recom- 
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bination of e + Fe^+ Fe+ is very small relative to the orbital 
timescale, the calculations should be valid. 

The recombination timescale is Tree = l/(arGc»T-e) seconds, 
where a^^ciT) is the recombination rate coefficient (in units of 
cm^ s~^) at temperature T and is the electron n umber density . 
For T ~ 10'' K, arec(T) ~ 4 x 10"^^ cm^ s'^ jNahailll997h . 
Since densities near the critical density of the line are the primary 
focus, Tie ~ 10^ cm~^; this gives Tree ~ 7 hours, which is much 
smaller than the orbital timescale, even during periastron passage, 
which takes approximately one month. Any possible time-delay ef- 
fects should not strongly affect the results as the light-travel-time 
to material at distances of ~ 0.5" in 77 Car is only about one week. 

3. 3. 1 Summary of the Basic Procedure 

Synthetic slit spectro-images are generated using a combination of 
Interactive Data Language (IDL) routines and radiative transf er cal- 
culat ions performed with a modified version of the SPLASH jPricd 
l2007 f) code. Here, we outline the basic procedure, with the specifics 
discussed in the following subsections. 

Nearly all of the required calculations are performed within 
SPLASH, which reads in the 3-D SPH code output for a specific 
orbital phase, rotates the data to a desired orbital orientation on the 
sky, and computes the line-of-sight velocity for all of the material. 
Next is the computation of the ionization volume created by r/B 
where Fe^ is photoionized to Fe^"*" . It is assumed that any material 
within this volume is in collisional ionization equilibrium, and that 
no material with T > 250, 000 K emits. All other material in the 
photoionization volume has its forbidden line emissivity calculated 
using Equation ( I13t , derived below. Material located outside the 
photoionization volume does not emit. 

The intensity / is computed by performing a line-of-sight in- 
tegration of the emissivity at each pixel, resulting in an image of 
the spatial distribution of the [Fe III] intensity projected on the sky. 
This is done for each velocity ii(,i„ used along the dispersion axis 
of the synthetic slit in the velocity range of interest, resulting in 
a set of intensity images. The individual [Felll] A4659 intensity 
images are then combined and convolved with the HST/STIS re- 
sponse using IDL routines in order to create a synthetic position 
versus velocity spectro-image for comparison to the observations. 

3.3.2 Defining the Binary Orientation Relative to the Observer 

With the two stars orbiting in the xy plane, the 3-D orientation of 
the binary relative to the observer is defined by the orbital inclina- 
tion i, the prograde direction angle 6 that the orbital plane projec- 
tion of the observer's line-of-sight makes with the apastron side of 
the semi-major axifl and the position angle on the sky of the +2: 
orbital axis, PA^ (Figure[3](. 

Following the standard definition for binary systems, i is the 
angle that the line-of-sight makes with the +2 axis. An i = 0° is 
a face-on orbit with the observer's line-of-sight along the +2 axis 
and the two stars orbiting counter-clockwise on the sky. An i = 90° 
places the line-of-sight in the orbital xy plane, while i = 180° is a 
face-on orbit with the two stars orbiting clockwise on the sky. 

For i = 90°, a value of = 0° has the observer looking along 
the +a; axis (along the semi-major axis on the apastron side of the 
system), while a.6 = 90° has the observer looking along the +y 

^ is the same as the angle (f> defined in figure 3 of O08. 




Figure 3. Diagrams illustrating the observer's position. Top: Schematic 
defining the inclination angle i that the observer's line-of-sight makes with 
the +z orbital axis, and the equatorial projection angle d of the line-of- 
sight relative to the apastron side of the semi-major axis x. The background 
orthogonal planes show slices of density from the 3-D SPH simulation at 
apastron in the xy orbital plane and the yz plane perpendicular to the or- 
bital plane and major axis. Bottom: Diagram defining the position angle on 
the sky, PAz, of the +z orbital axis (blue), measured in degrees counter- 
clockwise of north (N). The binary orbit projected on the sky is shown in 
yellow, as are the semi-major (red) and semi-minor (green) axes. The small 
inset to the left is a HST WFPC2 image of the Homunculus (Credit: NASA, 
ESA, and the Hubble SM4 ERO Team) included for reference. 



axis. In the conventional notation of binary orbits, the 'argument of 
periapsis' uj = 270° — 6. 

PAz defines the position angle on the sky of the +2 orbital axis 
and is measured in degrees counter-clockwise of N. A PAz =312° 
aligns the projected orbital axis with the Homunculus polar axis 
(Davidson etal. 2001; Smith 2006b, 2009) and has +2 pointing 
NW on the sky. A PA^ = 42° places the orbital axis perpen- 
dicular to the Homunculus polar axis with -\-z pointing NE. A 
PAz = 132° (222°) is also aligned with (perpendicular to) the 
Homunculus polar axis and has +z pointing SE (SW) on the sky. 
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3.3.3 Volume of Material Photoionized by rjB 

Based on the models of 77a's envelope by HOI and H06, it is as- 
sumed that Fe+ is initially the dominant ionization state of iron 
in primary wind material located r > 50 AU from tja, and that 
hydrogen and helium are initially neutral in the primary wind for 
distances of r > 155 AU and > 3.7 AU, respectively. The cal- 
culation of the volume of materi al photoionized from F e "^ to F e^^ 
by 77b follows that presented in iNussbaumer & Vogell ( 1 19871) for 
symbiotic star systems. For cases where hydrogen or helium alone 
determine the ionization structure, the boundary between neutral 
hydrogen (H'') and ionized hydrogen (H^) can be expressed ana- 
lytically. As r^A's extended wind consists mostly of hydrogen and 
helium (figure 9 of HOI), and because the [Felll] emission arises 
in regions where hydrogen is ionized but helium is not, such a cal- 
culation should provide a reasonable approximation for the size 
and shape of the photoionization boundary between Fe"*" and Fe^^. 
This, however, is an upper limit for the photoionization volume 
since the IP of 16.2 eV required for Fe^+ is slightly larger than 
the 13.6 eV needed to ionize H'\ 

In our calculation, the primary star is separated by a dis- 
tance p from a hot companion star that emits spherically sym- 
metrically Ln photons s^^ that are capable of ionizing H". The 
primary has a standard spherically symmetric mass-loss rate of 
Ml = 4:-Kr'^ jj.mu'n{r)voo, where n(r) is the hydrogen number 
density, jj, is the mean molecular weight, mu is the mass of a hy- 
drogen atom, and Voo is the terminal speed of the primary wind. 

In a small angle A9 around the direction 6, the equilibrium 
condition for recombination-ionization balance is 

AO fe 

Lh = A9 n(s)nJs)aB(}i,T^) ds, (1) 

471- Jo 

where s measures the distance from the ionizing secondary star 
in the direction 6, se is the boundary between H" and H"'', as 
is the total hydrogenic recombination coefficient in case B at 
electron temperature Te, an d Ue is the electron number density 
jNussbaumer & Vogellll987l) . The electron density can be written 
as 



In the modified SPLASH routine. Equation l[4j is solved for a 
specified Lh . However, since Equation l[4j does not account for the 
effects of the WWC, we assume that mass is conserved and that all 
primary wind material normally within the cavity created by rjB is 
compressed into the walls of the WWC region, with the ionization 
state of the WWC zone walls before being photoionized the same 
as that of the primary wind (i.e. iron is in the Fe"'" state). 

The stellar separation is computed directly from the 3-D SPH 
simulation. We use the same value for the abundance by number 
of He to H as HOI and H06, a(He) = 0.2. A constant value of 
QB = 2 . 56 X 10 "^^ cm^ at Te = 10,00 K is also used 
(Osterbrock 1989). Based on lMehner et alJ l l2O10l) and lVemer et al.l 
(2005), the value of Lh is chosen to be that of an 05 giant with 
Teff « 40, 000 K, which according to iMartins et al.l J2005h has 
a hydrogen ionizing flux of IQ*^ ** photons s~^. Any shielding 
effects due to the presence of the Weigelt blobs or other dense, 
slow-moving equatorial ejecta are not presently included. 

The third column of Figure |6] illustrates the time variability 
of the photoionization region created by tjb for slices in the xy 
orbital plane. The region is largest around apastron, when t^b is 
farthest from tja's dense wind. As rjB moves closer to t^a, thepho- 
toionization zone becomes smaller and more wedge-shapeqj due 
to the gradual embedding of tjb in 77a's wind. During periastron, 
rjB becomes completely enshrouded in 77a's thick wind, preventing 
material at large distances (> 10a) from being photoionized. As 
rjB emerges after periastron, the photoionization region is restored 
and grows as rjB moves back toward apastron. This plunging into, 
and withdrawal from, tja's wind by rjB leads to the illumination of 
distant material in very specific directions as a function of phase. 

3.3.4 The Emissivity Equation 

For a simple two-level atom, the volume emissivity j (erg cm~^ 
s^^ sr^^) of a forbidden line is 

j = -^hiyziNiAii , (7) 
47r 



ne(r) = (l-ha(He))n(r) , (2) 

where 

iirr^ firnHVoo 

and a(He) is the abundance by number of He relative to H. 
Expressing both s and r in units of p, defining u = s/p, 
and making the further approximation that qb is constant 
dNussbaumer & Vogellll987h . Equation ([T} can be written as 



where N2 is the number density of atoms in t he excited upper level 
and U 21 is the frequency of the line transition jDopita & Sutherland! 
l2003h . Using the total number density of element of interest E i n 
ionization state i, rii^E ~ Ni + N2 jlgnace &Brimeveill200^ . 
together with Equation (B5), one finds 

N2^( )ir^.,E~N2). (8) 

\n^q2i + A21 } 

Solving for N2 and using Equation (B4) for q2\ /qi2 gives 



where 



and 



X"^* = f{u, 6) , 



H+ Aixji^mB 



X" = 



r / °° 



aB(l + a{ile)y 



(a;2 - 2a;cos6i + 1)- 



dx . 



(4) 



(5) 



(6) 



91 



f hv2i \ 



N2 = rii^E <i 1 + — exp , „ , 
92 \kT J 



1 + ' 



(9) 



where ric, defined in Equation (B6), is the critical density of the 
line. By defining Qi,E = rn^E/nE as the fraction of element E 
in ionization state i, Ae = ue/un as the abundance of element 



* The outer edge looks circular only because this marks the edge of the 
spherical computational domain of the SPH simulation. 
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E relative to all nucl eons un, and 7e = un j n^. as the ratio of 
nucleons to electrons ilgnace & Brimeve3l2006h , rii^E becomes 



temperatures. The model therefore assumes zero [Felll] emission 
from material with T > 250, 000 K. 



ni,E = Qi.sAE'yenfi 

Therefore, using Equations ^ and J 10b . 



(10) 



47r 



1+^ 



(11) 



Based on Equation Jl lb . emission from a specific forbidden 
line is concentrated only in regions that (1) contain the right ion- 
ization state of the element of interest (in our case, Fe^"'"), (2) are 
near the critical density of the line (for [Fe III], ~ 10^ cm~^), and 
(3) are near the right temperature (for [Felll], T fa 32, 000 K). 
This important point is crucial to understanding the forbidden line 
emission observed in 77 Car. 

To compute the emissivity using SPLASH, the line profiles are 
assumed to have a Gaussian thermal broadening, so that Equation 
dUb is weighted by an exponential with the form 



hV2lA2lQi.EAEfene 

J = -, exp 



X < 1 + exp 



Vbi' 



Vlos 
Vth 



hV21 

kT 



31 ^ AaipiT^ 

32 fiQ,12ne 



, (12) 



where Equation (B8) has been used for ric, vios is the line-of-sight 
velocity of material in the 3-D SPH simulation for a specified ori- 
entation of the binary orbit relative to the observer, vt^ is a thermal 
velocity dispersion (taken to be 25 km s~^), and Vbin is the bin size 
used along the dispersion axis of the synthetic slit in the wavelength 
range of interest, expressed in thermal velocity units (25 km s~^). 

To relate the physical density in the SPH simulation to the 
electron number density, we use = Psph/ifJ'emn), with fie the 
mean molecular weight per free electron. As the focus is on forbid- 
den lines of trace metals, the value of rie should not be impacte d 
by the ionization balance of the metals ( Ignace & Brimeve3l2006h . 
Moreover, since far-UV radiation (or high-energy collisions) is nec- 
essary for the formation of the high-ionization forbidden lines, 
wherever [Fe III] emission occurs, hydrogen should also be ionized. 
For such regions dominated by H"*", we thus take /ie = 7e = 1. The 
volume emissivity then takes the final general form. 



J = 



hU2lA2lQi,EAEp. 



ph 



exp 



Vbi' 



Vlos 
Vth 



1 + exp 



hi'21 
kT 



gi ^ A2igiT ' mn 
32 



l3Ql2psph 



.(13) 



3.3.5 The Intensity and Synthetic Slit Spectro-Images 

Since the [Fe III] emission is optically thin, the intensity / is simply 
the integral of the volume emissivity j along the line-of-sight. 



jdl 



(14) 



( lMihalasll978h . The intensity is computed in SPLASH by perform- 
ing a line-of-sight integration of j through the entire 3-D simulation 
at each pixel, resulting in an image of the intensity in the [Fe III] 
A4659 line centered at a particular velocity (wavelength), as an ob- 
server would see it projected on the sky. This is done for multiple 
velocities {vbin) in 25 km intervals for line-of-sight velocities 
from -600 km s"^ to +600 km 

The resulting series of [Fe III] intensity images are combined 
using IDL routines to create a synthetic position versus velocity 
spectro-image. The IDL code reads in the SPLASH images and ro- 
tates them to a specified PA on the sky corresponding to a desired 
PA of the HST/STIS slit, assuming that the orbital axis is either 
aligned with the Homunculus polar axis at a PA2 = 312°, or ro- 
tated relative to the Homunculus axis by some specified angle. Each 
image is cropped to match the 0.1" width of the STIS slit, with the 
slit assumed centered exactly on the 77 Car central source. Since po- 
sition information is only available along one direction, each image 
has its intensity values integrated along each row of pixels within 
the slit width. This produces a single 'slice' of the [Felll] A4659 
intensity along the slit, centered at a specific velocity. These slices 
are combined to create the synthetic spectro-image. 

The resulting model spectro-images have spatial 0.003") 
and spectral (25 km s~^) resolutions that are better than the 
HST observations. The synthetic spectro-images are therefore con- 
volved with the response of HST/STIS in order to match its spa- 
tial (0.1") and spectral (37.5 km s^^) resolutions. Point-spread 
functions (PSFs) for STIS generated using the Tiny Tim program 
l lKrist&Hoo"3l 19991) are used for the convolution in the spatial di- 
rection, while a gaussian is used for the convolution in the spectral 
direction. Color in the model spectro-images scales as the square 
root of the intensity, with the colorbar ranging from to 1/4 of the 
maximum intensity, the same as used for displaying the observa- 
tions. 

Figure |4] shows that there is a noticeable difference be- 
tween the unconvolved (left) and convolved (right) model spectro- 
images, especially in the central ±0.15". Details in the uncon- 
volved spectro-images that arise near the inner WWC zone are 
completely unresolved in the convolved images, resulting in a 
bright, central streak of emission extending from negative to pos- 
itive velocities. This is consistent with the observations and indi- 
cates that HST/STIS lacks the spatial resolution needed to resolve 
the details of rj Car's inner WWC zone. 



The appropriate value of Qi^E is found assuming collisional 
ionization equilibrium, with the fr action of Fe^^ as a f unction of T 
based on the ion fraction data from lBrvans et al.l j2006l) . These data 
show that at T ~ 250, 000 K, the fractional abundance of Fe^"*" 

IS 

10-5-453^ with the remaining Fe being collisionally ionized to 
higher states. The amount of Fe^^ decreases even more at higher 



4 RESULTS FROM THE 3-D DYNAMICAL MODEL 

In the discussions below, it is assumed for simplicity that phase zero 
of the spectroscopic cycle (from the observations) coincides with 
phase zero of the orbital cycle (periastron passage). In a highly- 
eccentric binary system like 77 Car, the t wo values are not e xpected 
to be shifted by more than a few weeks jGroh et alj20I0ah . Such a 
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Unconvolved 



Convolved with HST Response 




Figure 4. Example synthetic spectro-images of [Fe III] A4659 both uncon- 
volved (left) and convolved with the response of HST/SJIS (right), assum- 
ing i = 138°, e = 0°, and PA^ = 312°. Roman numerals I and II 
indicate the top and bottom of the slit, respectively. Color is proportional to 
the square root of the intensity, and the velocity scale is ±600 km s~^. 



time shift would only cause a small change of ~ 10° in the derived 
best value of 9, which will not affect the overall conclusions. 



4.1 Hydrodynamics of the Extended WWC Region 

The first two columns of Figures [5] and [6] show, respectively, the 
density and temperature in the xy orbital plane from the 3-D SPH 
simulation. The binary system has undergone multiple orbits, as in- 
dicated by the dense arcs of wind material in the outer (> 20a) re- 
gions. On the ~x (periastron) side of the system are narrow cavities 
carved by tjb in tja's dense wind during each periastron passage. 
The cavities are most easily seen in the temperature plots as they 
contain warm (~ lO** — 10^ K), low-density wind material from 
77b. Bordering these cavities are compressed, high-density shells of 
primary wind that form as a result of the WWC. 

The insets of Figure [6] illustrate the creation of one of these 
narrow wind cavities and its bordering dense shell of primary wind. 
Their spiral shape is due to the increased orbital speeds of the stars 
during periastron. This is in contrast to phases around apastron 
(rows b, c and d) when orbital speeds are much lower and the cur- 
rent WWC zone maintains a simple, axisymmetric conical shape 
(O08; P09; Pll). The increasing orbital speeds when approaching 
periastron causes the postshock gas in the leading arm of the WWC 
zone to be heated to higher temperatures than the gas in the trailing 
arm (insets of rows d and e), an affect also found by PI 1. 

Following periastron, rjB moves back to the +x (apastron) 
side of the system, its wind colliding with and heating the dense 
primary wind that flows unimpeded in the +x, —y direction. As 
pointed out by PI 1, the arms of the WWC region shortly after this 
{(j) « 0.1) are so distorted by orbital motion that the leading arm 
collides with the trailing arm from before periastron, leading to ad- 
ditional heating of the postshock gas in the trailing arm (compare 
cf> ~ 1.122, Figure |6] row f). The leading arm of the WWC zone 
(including the portion that collides with the trailing arm) helps to 
form a dense, compressed shell of primary wind that flows in the 
-\-x and — y directions after periastron passage. Because the wind of 
77b collides with a dense wall of postshock primary wind with high 
inertia, the primary wind controls the overall rate of expansion of 
the resulting spiral (Pll). 

The overall stability of the expanding shell of primary wind 
on the apastron side of the system depends on the shock thickness 
(Vishniac 1983; Wunsch et al. 2010; Pll). Portions of the shell 
moving in the — j/ direction appear to be the most stable due to the 



increased amount of primary wind that borders it in this direction. 
However, our simulations show that at (ji ~ 0.3, the upper portion 
of the shell expanding in the +x, +y direction starts to fragment. 
Eventually, the wind of rjB is able to plough through the shell, caus- 
ing it to separate from the leading arm of the WWC zone (</!> f» 0.4, 
row c of Figure|6]l. This produces a pair of dense 'arcs' of primary 
wind on the apastron side of the system. The arc expanding in the 
+x, —y direction is the remnant of the leading arm of the WWC 
region from the previous periastron passage. The arc in the +x, +y 
direction is the remnant of the trailing arm of the WWC region from 
just prior to the previous periastron passage. Multiple pairs of these 
arcs are visible in Figures[5]and[6l They are also quite spatially ex- 
tended; by the time the system is back at periastron, the arcs from 
the previous periastron are up to 70a from the central stars (in the 
xy plane). As the arcs expand further, they gradually mix with the 
surrounding low-density wind material from 77B . 

The fragmentation of the shell and formation of the arcs is 
mainly due to the above-mentioned collision of the leading arm of 
the WWC zone with the trailing arm just after periastron. This col- 
lision produces instabilities at the interface between the two arms 
(Pll), which, together with the shell's expansion, the pressure from 
t/b's high- velocity wind, and the lack of bordering primary wind 
material for support, causes the shell to break apart in the +x, +77 
direction ~ 1.7 years after periastron. 

While the term 'shell' is used to describe the layer of com- 
pressed primary wind that flows in the direction of apastron fol- 
lowing periastron passage, this is not related t o the 'shell ejection' 
event discussed in the past rj Car l i terature 1 Zanella et alj 19841; 



Davidsojl999l ; lDavidson et al.ll999l;ISmith et al.[2003d ; lDavidsonl 
200I). T he scenario prop osed by Izanella et alj ^19840 and advo- 
cated bv' Pavidsoij ([l999) to explain 77 Car's spectroscopic events 
is a qualitative single-star model wherein some sort of thermal or 
surface instability is presumed to induce significant mass loss from 
7; Car approximately every 5.54 years. A shell ejection as the result 
of a latitude-dependen t disturbance in the wind of a sin g le star has 
also been proposed ( Smith et alj2003j ; [Dav idson 20051) . iDavidsorj 
( 2005) suggested that the shell event may be triggered by the close 
approach of a secondary star, but this is still mainly a single-star 
scenario requiring some kind of surface instability in the primary. 

In contrast to these ideas, the formation of the outflowing shell 
of primary wind seen in the 3-D SPH simulations is a natural, un- 
avoidable consequence in any high-eccentricity, massive CWB in 
which the primary star has a significant mass-loss rate (significant 
both in terms of overall mass-loss rate and mass-loss rate relative 
to the companion star). It is not so much an 'ejection' event as it 
is a chance for the normal primary wind to flow in the direction 
of apastron for a brief time (~ 3 — 5 months) while 773 performs 
periastron passage. Eventually, 773 returns to the apastron side and 
its high-velocity wind collides with and compresses this primary 
wind material, forming a dense shell that continues to propagate in 
the direction of apastron until it breaks apart at its weakest point, 
resulting in a pair of dense 'arcs', as described above. 

4.2 Physical Origin and Location of the Broad 
High-Ionization Forbidden Line Emission 

Figure |5]presents slices in the orbital plane from the 3-D SPH sim- 
ulation at (/> = 0.976, used to model the blue-shifted emission arcs 
seen at STIS slit PA = +38°. The colors show log density, log 
temperature, photoionization zone created by r7B, and square root 
of the modeled emissivity of the [Fe III] A4659 line. The photoion- 
ization region is confined to the same side of the system as 77B , the 
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Figure 5. Snapshots in the orbital xy plane from the 3-D SPH simulation of r) Car aX (f) = 0.976 used to model the observed blue-shifted emission arcs of 
Figure[2] Color shows, from left to right, log density, log temperature, photoionization volume created by rj^ (white = ionized), and square root of the modeled 
emissivity of the [Fe III] A4659 line. The color bar in the last panel has been adjusted to make faint emission more visible. The box size is ±105a ±1622 AU 
~ ±0.7" {D = 2.3 kpc). Axis tick marks correspond to an increment of 10a 154 AU ~ 0.067". In the last panel, [Fe III] emission only originates from 
material within the photoionization volume that is near the critical density of the [Fe III] line (~ 10^ cm""^) and at the appropriate temperature. 



dense wind of t^a preventing the far (—a;) side of the system from 
being photoionized. All of the [Fe III] emission is thus concentrated 
in two spatially-distinct regions on the apastron side of the system: 

• The strongest [Fe III] emission originates near the walls of the 
current WWC zone in the inner ~ 30a « 0.2" of 77 Car. 

• Faint, spatially-extended (out to ^ 60a ~ 0.4") [Felll] emis- 
sion arises in the arcs of dense primary wind formed during the 
previous periastron passage. 

The [Fe ill] emission occurs in these two areas because they are the 
regions within the photoionization zone that are near the critical 
density of the [Fe III] line and at the appropriate temperature. 

Figure [6] illustrates how the density, temperature, photoion- 
ization region, and [Felll] A4659 emissivity change with orbital 
phase. The six phases shown correspond to those observed in figure 
11 of G09. During periastron passage {(j> — 0.045 and 1.040, rows 
a and e), i]b becomes deeply embedded in the dense wind of t^a 
and the photoionization region shrinks considerably, down to only 
a few semi-major axes across. This leads to an effective 'shutting- 
off' of the forbidden line emission at large distances (> 10a). By 
(j) = 0.213and 1.122 (rows b and f), rye has moved far enough from 
rjA to re-establish the large photoionization zone and spatially- 
extended forbidden line emission. The photoionization region is 
largest, and the [Felll] most spatially extended, around apastron 
(row c, (f> — 0.407). By = 0.952 (row d), r^B starts to become 
embedded in the dense wind of r^A, the photoionization region in 
the +x, —y quadrant decreasing slightly in size. 

Far-UV radiation from rjB leads to highly ionized regions that 
extend outward from its low density wind cavity in the direction of 
system apastron. During periastron passage, the disappearance of 
the extended high-ionization forbidden emission can be attributed 
to the wrapping of the dense primary wind around rye, which traps 
its far-UV radiation and prevents it from photoionizing the outer 
wind structures responsible for the observed emission. While tjb is 
embedded in 77a's wind, the latter flows unimpeded in the direction 
of apastron, eventually forming dense arcs of material that expand 
outward as t^b completes periastron passage and the inner regions 
of the WWC zone regain their near conical shape. These cooler, 
dense arcs of expanding primary wind also produce high-ionization 
forbidden line emission when photoionized by rjB- 



Observed temporal variations of the high-ionization forbidden 
lines can thus be linked to the orbital motion of rjB in its highly 
eccentric orbit, which causes different portions of the WWC re- 
gions and extended arcs of primary wind from earlier cycles to 
be photoionized. As the two stars move closer to or farther from 
each other, orbital motion also leads to changes in the density and 
temperature of the inner WWC zone, which in turn modifies the 
size and shape of ?7b's photoionization volume, and thus the over- 
all shape, location, and intensity of the high-ionization forbidden 
line emission, even at phases away from periastron (rows b, c, d, 
and f of Figure [6ll. Forbidden line emission from the compressed, 
inner WWC region increases in intensity until the critical density 
is reached or the temperature exceeds that at which the appropriate 
ions can exist. Forbidden emission from the extended arcs of pri- 
mary wind decreases in intensity as the arcs gradually expand and 
mix with the surrounding low-density wind material from ijB ■ 



4.3 Synthetic Slit Spectro-images and Constraining the 
Orbital Orientation 



4.3.1 Constraint 1: Emission Arcs at Slit PA = +38° 



0.976 



The blue-shifted emission arcs in Figure[2]represent distinct, well- 
defined structures observed at a specific orbital phase and slit PA. 
As such, they provide a natural basis for modeling and can be used 
to constrain the orbital orientation. We have performed a parameter 
study in i, 6, and PA^ with the goal of determining which set(s), 
if any, of orientation parameters result in synthetic spectro-images 
that closely match the observations in Figure[2] The value of d was 
varied in 15° increments for values of 0° ^ S ^ 360°, with i varied 
in 5° increments over the range 0° ^ i ^ 180° for each 9 value, 
and PAz varied in 5° increments over the range 0° ^ PAz ^ 360° 
for each pair of i, 9 values. 

We find that only synthetic spectro-images generated for 
— 15° ^ S ^ +45° are able to reasonably match the observations 
taken al <j) = 0.976, PA = 38°. However, an ambiguity exists in 
i and PA^, with values of 30° ^ i < 50°, PA^ = 272° to 332° 
(top row of Figure|71l and 130° sC i ^ 150°, PA^ = 282° to 342° 
(bottom row of Figure [T} both capable of producing entirely blue- 
shifted arcs that resemble those observed. Our best morphological 
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Figure 6. Same as Figure|5] but for = 0.045, 0.213, 0.407, 0.952, 1.040, and 1.122 (rows, top to bottom), which correspond to phases at which HST/STIS 
observations were accomplished (figure 1 1 of G09). Insets in the first two columns are a zoom of the inner 10a, included to illustrate the complex dynamics 
of the 'current' WWC region. 
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Figure 7. 2-D spatial distribution on the sky of the square root of the modeled intensity in the [Fe III] A4659 line and synthetic spectro-images for orbital 
orientations i = 42°, d = 7°, PA^ = 302° (top row) and i = 138°, 9 = 7°, PA^ = 317° (bottom row), which lie near the centers of the two derived 
best-fit ranges of orientation parameters for matching the observations taken at = 0.976, PA = 38° in Figure |2] Columns are, from left to right: 2-D 
distribution on the sky of the square root of the modeled intensity in the [Felll] A4659 line for blue-shifted material with vi^s between —500 km and 
— 100 km s~^; same as first column, but for red-shifted material with vi^s between +100 km and +500 km s~^; model spectro-image convolved with 
the response of HST/STIS; and observed spectro-image with mask. The projected x, y, and z axes are shown for reference in the 2-D images of the intensity 
on the sky, as is the direction of north. The 0.1" wide STIS slit at PA = +38° is also overlaid. Roman numerals I and II indicate the top and bottom of 
the slit, respectively. All following 2-D projections of the modeled [Fe III] intensity in this paper use this same labehng convention. The color scale in the 
spectro-images is proportional to the square root of the intensity and the velocity scale is ±600 km s~^. Lengths are in arcseconds. 



fits are obtained using values of i = 42°, e = 7°, PA^ = 302° 
or i = 138°, 9 = 7°, PAz = 317°, which lie near the centers of 
these two derived ranges of best-fit orientation parameters. 

The observed morphology and asymmetries in brightness and 
vio3 of the spatially-extended emission arcs determine the allowed 
ranges of the orbital orientation parameters. The parameters i and 
PAz primarily control the brightness asymmetry and spatial orien- 
tation of the blue- and red-shifted emission components on the sky. 
Figure|7]shows that the synthetic spectro-images at slit PA = 38° 
for i = 42° and 138° (at nearly identical 9 and PAz) are remark- 
ably similar. Yet, the two orbital orientations are drastically dif- 
ferent, resulting in unique distributions on the sky of the blue- and 
red-shifted emission. The red component is the most notable, which 
extends to the SE for i = 42°, but to the NW for i = 138°. The 
asymmetry in brightness between the two blue-shifted arcs is also 
controlled by i. Only values of 30° ^ i ^ 50° or 130° sC i < 150° 
produce a brightness asymmetry, with other i resulting in arcs that 
are approximately equal in brightness. 

The value of PAz determines the final orientation of the pro- 
jected z axis, and thus of the individual emission components. Be- 
cause the lower blue-shifted arc in the observations is brighter than 
the upper arc, the extended blue emission projected on the sky must 
be brighter in directions to the SW (where the lower half of the slit 
is located, position II in Figure[7}, and dimmer in directions to the 
NE (position I). This limits PAz to values between 272° and 342°. 

IVIore importantly, i and PAz determine the 3-D orientation 
of the orbital axis. If i = 42° and PAz = 272° to 342°, the 
orbital axis is inclined toward the observer and is not aligned 
with the Homunculus polar axis. The two stars would also orbit 
counter-clockwise on the sky. However, if i = 138°, then, because 



PAz ~ 312°, the orbital axis is closely aligned in 3-D with the 
Homunculus polar axis, and the stars orbit clockwise on the sky. 

Figure [8] shows that the value of 6 (for fixed i = 138° 
and PAz = 312°) determines the vios of the material exhibit- 
ing [Felll] emission. According to our 3-D model, the spatially- 
extended [Fe III] emission originates in the expanding arcs of dense 
primary wind formed during the previous periastron passage. When 
9 ~ Q° (top row of Figure [8]l this emission is mostly blue-shifted 
because the emitting material within t/b's photoionization zone 
is moving mostly toward the observer. The model spectro-image 
consists of spatially-extended, blue-shifted arcs because the slit at 
PA = 38° primarily samples the blue component, which stretches 
from NE to SW on the sky. In contrast, the extended red component 
is not sampled; only a small amount of red-shifted emission in the 
very central ±0.1" core falls within the slit. Moreover, only values 
of —15° 55 ^ ^ +45° produce significant amounts of spatially- 
extended (~ 0.7" in total length) blue-shifted emission with the 
correct observed asymmetry in vios, wherein the dimmer upper arc 
extends ~ 75 km s^^ more to the blue than the brighter lower arc. 

Figure|7]demonstrates that our 3-D dynamical model and de- 
rived best-fit orbital orientation(s) reproduce all of the key features 
in the observations taken at = 0.976, PA = 38°. Both syn- 
thetic and observed spectro-images contain concentrated, velocity- 
extended emission in the central ±0.1" core. The completely blue- 
shifted emission arcs extend to roughly the same spatial distances 
(~ ±0.35"). The asymmetry in brightness between the upper and 
lower arcs is matched as well. The asymmetry in vios is also repro- 
duced and of the same magnitude. However, the arcs in the model 
images stretch a bit farther to the blue, to ~ —550 km (upper 
arc) and —475 km s^^ (lower arc), versus ~ —475 km and 
—400 km s^^ in the observations. The synthetic images thus show 
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Figure 8. Same as Figure|7] but for values of 6 = 0° , 90° , 180° , and 270° (rows, top to bottom), assuming i = 138° and PAz = 312° . 



emission at velocities slightly above the value used in the SPH sim- 
ulation for j7a's wind terminal speed. This is a minor discrepancy 
though, likely due to the extended emitting material receiving an 
extra 'push' from the fast wind of rjB (see subsection 4.3.4). 

4.3.1.1 Results for Orbital Orientations with a f» 180° 
Returning back to Figure [S] one sees that the spectro-image for 
6 = 180° fails to match the observations. There is a distinct lack 
of any spatially-extended, blue-shifted emission. Instead, the ex- 
tended emission arcs are entirely red-shifted due to the quite dif- 
ferent vios of the emitting material. When 9 = 180°, riA is be- 
tween the observer and r]B during most of the orbit. Therefore, all 
of the emitting material within tjb 's photoionization zone is on the 
far side of the system and moving away from the observer. The ex- 
tended forbidden line emission is thus mostly red-shifted. The slit at 
PA = 38° now primarily samples this red component, producing a 
spectro-image consisting of entirely red-shifted, spatially-extended 
arcs, the direct opposite of the observations. 

Changing the value of i and/or PAz when 9 — 180° does 



not result in a better match to the observations since most of 
the material photoionized by rjB and exhibiting forbidden line 
emission is still moving away from the observer. This is true 
for all values of 6 near 180°; none of the model images gener- 
ated for 135° ^ S ^ 225°, regardless of the assumed i and 
FAz, matches the (f> = 0.976, PA = 38° observations. Ex- 
ample model spectro-images for these orientations can be found 
in [Madural jzOlO ). We therefore find that an orbital orientation 
that places rjB on the ne ar side of ?7a at pe ri astron, such as 



that fav o red by Falceta-G oncalves et alj (1200: 



Abraham et al 



(2005b); Abraham & Falceta-Goncalves ( 2007); Kashi & Sokei 



(2007, 2008); Falceta -Goncalves & Abraham C2009.) and others, is 
explicitly ruled out. 

4.3.1.2 Results for Orientations with 6 ^ 90° and 270° 

Model spectro-images for 6 = 90° and 270° also fail to match 
the observations. The 2-D images of the modeled [Fe III] intensity 
on the sky in Figure[8]show that in both cases, one half of the STIS 
slit is empty (from center to position II for 9 = 90° and from I 
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Figure 9. Same as Figure|2] but for = 0.952, slit PA = —28°. Note that the red- and blue-shifted emission in the model spectro-image that assumes 
i = 42° (top row) spatially extends in the wrong directions compared to the observations, r) Car's orbital orientation parameters are thus constrained to values 
of j Si 130° to 145°, Si 0° to 15°, PA^ si 302° to 327° (see text). 



to the center for 6 = 270°). Therefore, emission is absent to one 
side spatially in the spectro-images. Changing the value of i has no 
effect on this since a change in i at these 6 values is equivalent to 
a rotation about the x axis in the 2-D intensity images. Instead, a 
rotation of nearly 90° in PAz is needed in order to place emitting 
material in both halves of the slit. However, such a rotation does 
not result in a match to the observations. Synthetic spectro-images 
for 45° < e < 135° and 225° < 6 < 345° all suffer from the 
same problems, namely, either no blue-shifted ring-like emission 
feature (if PAz ~ 132° or 312°), or ar c-like emissio n to only one 
side spatially (if PA^ 42° or 222°) jMadurall201ol) . This is true 
regardless of the i and PAz assumed. 



4.3.2 Constraint 2: Variations with Phase at Slit PA = —28° 

Using the 3-D dynamical model we generated synthetic spectro- 
images for a slit PA — —28° at each of the phases in figure 11 
of G09, </> = 0.045, 0.213, 0.407, 0.952, 1.040, and 1.122. The 
main goals were to determine if the binary scenario and 3-D model 
could explain the phase dependence observed in the high-ionization 
forbidden lines, and whether the set of PA — —28° observations 
could resolve the ambiguity between the sets of i = 35° — 50° and 
i — 130° — 145° best-fit orientations found in the above subsection. 

Figure |9] illustrates the results of investigating the two best- 
fit i regimes. It is clear that synthetic spectro-images for i fa 30° 
to 50° do not match the observations taken at slit PA = —28° 
(top row of Figure |9j. There are a number of discrepancies, but 
the most important and obvious is that the blue- and red-shifted 
emission components in the synthetic spectro-image extend in the 
wrong spatial directions, with the blue emission stretching to the 
NW and the red to the SE, the opposite of the observations. 

In contrast, synthetic images generated for i ~ 130° to 150° 
(bottom row of Figure|9) are capable of matching the observations. 
The morphology of the PA — —28° observations breaks the de- 
generacy in i and fully constrains, in 3-D, r] Car's orbit. Because 
red-shifted emission is observed to extend spatially in directions to 



the NW, the red emission component of [Fe III] on the sky must be 
oriented such that it too stretches NW and falls within the top half 
of the PA — —28° slit. The synthetic images of the [Felll] inten- 
sity on the sky (first two columns of Figures |9] and [TOt illustrate 
how the slit at PA = —28° samples the emitting regions of rj Car's 
extended interacting winds in very different directions compared to 
slit PA = 38° . This is why the observed spectro-images are so dif- 
ferent between the two slit PAs, even at similar orbital phases. We 
find that only orbital orientations with i ~ 130° to 145° , S ~ 0° 
to 15°, PAz ~ 302° to 327° are able to simultaneously match the 
observations taken at both slit PA = 38° and —28°. 

Figure[TO]shows that the 3-D dynamical model and derived or- 
bital orientation reproduce the overall observed shape and velocity 
structure of the [Fe III] emission at each phase for slit PA — —28°. 
Fi gurefTOlf urther illustrates how rj Car's extended interacting winds 
change with time. The orientation of the blue component does not 
appear to change, always stretching from NE to SW on the sky. 
However, the spatial extent does change, growing larger going from 
periastron to apastron as the wind structures flow outward. Moving 
from apastron back to periastron, the intensity and spatial extent of 
the blue component eventually start to decrease as the expanding 
arcs of primary wind drop in density and mix with the surrounding 
wind from rjB. The red-shifted emission changes in a way similar 
to that of the blue component, but points mainly NW. 



4.3.3 Constraint 3: Observed Variations with (f) and Slit PA 

Synthetic spectro-images were generated at the ten combinations of 
phase and slit PA in Table |2l assuming the best-fit orientation i — 
138°, 6* = 7°, and PA^ = 317°. The goal was to determine how 
well the 3-D model and derived binary orientation could reproduce 
observations at a variety of other phases and slit PAs. Model images 
for (j) = 1.001 and 1.013 are nearly identical to those at(j> ~ 0.045 
and 1.040 in Figure[lO]and add no new information. We therefore 
focus this discussion on the phases before and after periastron. 

Figure[TT]presents the results. The overall match between syn- 
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Figure 10. Same as Figurel?] but for phases (rows, top to bottom) = 0.045, 0.213, 0.407, 0.952, 1.040, and 1.122, with the STIS slit at PA = -28° (as in 
figure 11 of G09), and assuming i = 138°, 6 = 7°, PA^ = 317°. 
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Table 2. Phases and slit PAs from figures 12 and 13 of G09 modeled in 
Section 4.3.3 and shown in Figure fTTI 






PA 


Figure 


0.601 


+22° 




0.738 


-82° 




0.820 


+69° 


EH: 


0.930 


-57° 


EB 


0.970 


+27° 


[IB 


0.984 


+62° 


[nj 


0.995 


+70° 


[nt 


1.001 


+69° 




1.013 


+105° 




1.068 


-142° 





thedc and observed spectro-images is quite good. The model repro- 
duces all of the key spatial features, namely, extended (> ±0.1") 
[Fe III] emission that is almost entirely blue-shifted for positive slit 
PAs, but that is partially red-shifted for negative slit PAs. Observed 
variations in emission with phase are also reproduced. 

Model images at <;/> = 0.601 i fTTl i) and 0.738 (IB') demon- 
strate how a large change in slit PA results in very different spectro- 
images. The spatial distributions of the emission on the sky are very 
similar at these two phases, as expected since the system is near 
apastron when orbital velocities are their lowest. However, the slit 
at PA = —82° samples the emission in very different directions 
compared to PA = 22° , producing a spectro-image that resembles 
those taken at PA = -28° in FigurefTUl 

The partial ring of emission observed at = 0.984, PA = 
62° dl IF), is also noteworthy. The match between model and ob- 
servations is very good, the emission arc to the SW present in both, 
as well as the much shorter arc to the NE between ~ —500 km s^^ 
and —350 km s~^. Two effects appear to be causing the NE arc to 
be shorter. First, due to their outward expansion, the extended arcs 
of primary wind have dropped in density and started to mix with 
the surrounding low-density wind from ?7b. Second, the photoion- 
ization of material in directions to the NE by 77B has diminished 
due to its clockwise orbital motion and gradual embedding in »7a's 
wind. These lead to a decrease in the amount of extended [Fe III] 
emission to the NE. This, combined with the slit PA, produces a 
partial arc of blue-shifted emission to the NE in the spectro-image. 

By = 0.995 (life), the extended emission has vanished, im- 
plying that the ionizing flux of photons from t^b shuts off sometime 
between (j) = 0.984 and 0.995. At = 1.068 (Ilh), there is still no 
observed spatially-extended emission. Yet, the 2-D model images 
of the intensity on the sky show that the extended [Fe III] emission 
has started to return in directions to the NE. This is visible in the 
model spectro-image as a bright, blue-shifted bulge that points NE. 
Unfortunately, the observational data in the central ±0.15" is of 
insufficient quality to tell if such an emission bulge was detected. 
Nevertheless, the 3-D model predicts that tjb should start to emerge 
from ?7a's dense wind and begin to restore the spatially-extended, 
high-ionization forbidden line emission at </> ~ 1.068. 

4.3.4 Discrepancies Between the Observations and 3-D Model 

During most of the orbit, the wind of 773 collides with the dense 
arcs of primary wind on the apastron side of the system where 



the spatially-extended forbidden line emission forms (see rows b 
through d of Figure[6}. Because the arcs are bordered by a large, 
low-density wind cavity created earlier by t^b, they have almost no 
support against t^b's high- velocity wind, which is able to drive por- 
tions of the arcs into the cavity at velocities just above their outflow 
speed, the terminal speed of the primary wind. An orientation that 
places the observer on the apastron side of the system thus results 
in model spectro-images containing blue-shifted emission at veloc- 
ities slightly higher than the wind terminal velocity of the primary. 

Because the observations do not contain blue-shifted emission 
at speeds more than ~ 500 km s~^, the wind speed(s) of t^a and/or 
r/B used in the SPH simulation may be a bit too large. The launch- 
ing of the two winds at their terminal velocity in the simulation 
may also be having a small effect on the results. Better constraints 
on t/b's wind terminal speed and improved 3-D simulations with 
proper driving of the stellar winds should help resolve these is- 
sues. Improved modeling coupled with future observations of the 
extended forbidden emission in rj Car's central arcsecond also has 
the potential to further constrain the wind terminal velocity of 77A. 

The main discrepancies between synthetic spectro-images as- 
suming i = 138°, 9 = 7°, PAz = 317° and the observations 
for slit PA = —28° are: (1) the extended blue-shifted component 
is located ~ 100 km farther to the blue in the model images 
than in the observations, and (2) the red-shifted component in the 
model images is ~ 75 km broader than observed. Both issues 
are likely due to the reasons discussed above concerning the wind 
terminal velocities of rjA and -qs ■ 



5 DISCUSSION 

5.1 The Orientation and Direction of r] Car's Binary Orbit 

There has been much speculation about the binary orientation, and 
many papers published offering suggestions (see Section 1). The 
detailed modeling of Section 4 tightly constrains the observer's 
line-of-sight to angles of ^ ~ 0° — 15° prograde of the semi-major 
axis on the apastron side of the system, placing tjb behind rjA dur- 
ing periastron. Given the uncertainties in the stellar and wind pa- 
rameters of TjB used in the 3-D simulations, and the assumption that 
phase zero of the orbital and spectroscopic periods coincide, values 
of —15° < S < +30° may also be possible. A study in the stel- 
lar/wind parameters of tjb, or improved observational constraints, 
is needed to further refine the value of 6. 

More importantly, the results bound, for the first time, the 
orientation of rj Car's orbital plane, with the orbital axis closely 
aligned in 3-D with the inferred polar axis of the Homunculus at 
an i X 130° - 145° and PA^ ^ 302° - 327°. Figure [12] illus- 
trates the orientation of 77 Car's binary orbit on the sky relative to 
the Homunculus. With i f« 138° , 6' « 7° , and PA^ 317° , the re- 
sulting projected orbit on the sky has tjb moving clockwise relative 
to r/A, with TjB approaching r/A from the SW prior to periastron, 
and receding to the NE afterward. 

In their recent work modeling r) Car's RXTE light curve, O08 
and P09 investigated values of i in the range 0° < i < 90°, and 
each obtained a best-fit value of i ~ 42° , which they state is consis- 
tent with the orbital axis of the rj Car binary being aligned with the 
Homunculus polar axis. However, the inclination angle of the Ho- 
munculus, as defined in Davidson et al. (2001) and Smith (2006b|), 
corresponds to the tilt of the polar axis out of the line of sight (fig- 
ure 5 of Davidson et al. 2001). Since this is the same quantity as 
the inclination defined for binary orbits, the polar axis of the Ho- 
munculus is tilted from the plane of the sky by ~ 48° . Therefore, 
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Figure 11. A Same as Figure[To] but for the observed orbital phases and STIS sht PAs hsted in Table|2] Phases are (top to bottom): <j> = 0.601, 0.738, 0.820, 
and 0.930. The 0.1" wide STIS slit is shown overlaid at the appropriate PA in the 2-D images of the modeled intensity on the sky. Roman numerals I and II 
indicate the top and bottom of the slit, respectively. 



an i « 138° is required in order for the rj Car orbital axis to be 
aligned with the Homunculus polar axis, not 42°. 

The spatially-unresolved nature of the X-ray observations and 
overall top-bottom syrmnetry of the problem are what allowed O08 
and P09 to reasonably reproduce the RXTE light curve using an i ~ 
42°. However, an i « 138° works as well. Moreover, with either i, 
any PA on the sky of the orbital axis produces a good fit to the RXTE 
light curve. Unfortunately, this ambiguity in i and PA of the orbital 
axis has not been sufficiently discussed in the rj Car literature, and 
most have simply assumed that the orbital and Homunculus polar 
axes are aligned. Yet, knowing the value of i is crucial for proper 
interpretation of spatially-resolved observational diagnostics and, 
most importantly, determining the stellar masses. 

Taking into account this ambiguity in i inherent to 3-D mod- 
els of X-ray data, our derived ranges of i and 9 agree very well 
with those of O08, P09, an d PI 1. They are als o consistent with the 
i and uj values derived bv lGrohetalltolObl) in their attempts to 



understand the origin of the high- velocity (up to —1900 km s~^) 
absorption wing detected in He I A10830 during rj Car's 2009 
periastron passage. Our binary orientation is in complete agree- 
ment with simple models proposed to explain the available ra- 
dio (Du ncan & White 2003; White et al. 2005) and Hell A4686 
jTeodoro et al.l201 1\) observations as well. 

Further support for ou r results an d derived orbital orientation 
comes from recent work bv lMehner e t al. (2010), who mapped the 
flux of the broad, blue-shifted component of [Felll] A4659 and 
found that most of the emission originates in the inner 0.15" re- 
gion of rj Car (see their figure 8). The synthetic spectro-images and 
2-D spatial maps of [Fe III] A4659 in Section 4 show that most of 
the emission should indeed originate in the inner 0.15". According 
to our 3-D model, this emission forms in/near the current WWC 
zone in regions photoionized by rjB and manifests itself in spectro- 
images as a bright, central streak in the inner ±0.15" that spans a 
wide range of velocities from blue to red. This emission appears 
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Figure 11. B Same as Figure[TT]A, but at phases (top to bottom) = 0.970, 0.984, 0.995, and 1.068. 



as a bright central streak because HSTiSTIS lacks the resolution 
necessary to resolve the inner WWC region's complex features. 

lMehneretalJ ( l201(]|) additionally found that the spatial distri- 
bution of the blue-shifted component of [Fe III] A4659 in the cen- 
tral 0.1" of rj Car is not as sharp as a stellar point source and is 
detectably elongated to the NE and SW (see their figures 8 and 9). 
Model images displaying the 2-D spatial distribution on the sky of 
the blue-shifted component of [Fe III] A4659 (first column of Fig- 
ures □ [Tl and[TTJ show that the emission is elongated to the NE 
and SW for our suggested orbital orientation, even in the central 
0.1". While detailed mappi ng with HST /STIS is needed to verify 
the results of iMehner et al.l (12010), they do strongly suggest that 
our 3-D model and binary orientation are correct. 

All other orbital orientations, such as that with 6 ~ 180° fa- 
vored by Falceta-Goncalves et al. (2005); Abraham et al. (2 C)05b,) ; 
Abraham & Falceta-Goncalves ( 2007); Kas hi& Sokeri j2007l 
200&) ; . Falceta-Goncalves & Abraham, (i2009,) ; and iKashi & SokeJ 
( l2009h . produce spectro-images that are in strong disagreement 
with the (j) — 0.976, PA = 38° observations, lacking any 



spatially-extended, blue-shifted emission and containing signifi- 
cant amounts of extended, red-shifted emission. With t^b the main 
source of photons capable of producing high-ionization forbidden 
line emission, there is no realistic situation with 9 « 180°, for any 
i and PA^, that can produce spectro-images like those observed 
since there will always be an extended red-shifted component. 
Hypothesizing some method to absorb or diminish the red-shifted 
emission does not help since there would then be almost no 
spatially-extended emission detected at slit PA — 38°. 

There is no known source of ionizing photons in the 77 Car 
system that permits a ^ ~ 180° orientation. Photons generated in 
the WWC region will not work simply because the WWC zone is 
on the opposite side of the system during most of the orbit and t^a's 
dense wind would easily absorb any such photons before they reach 
the periastron side. Even if some contrived process could produce 
significant amounts of spatially-extended, blue-shifted emission on 
the periastron side of the system, this would not result in a match to 
the observations as the spectro-images would then contain both red- 
and blue-shifted extended emission. There would have to be some 
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Based on proper motions and observed blue-shifted emission, the 
Weigelt blobs are located on the same side of r/ Car as the ob- 
server and are thought to have been ejected sometime around the 
smaller eruption of 1890 th at also formed the Li ttle Homuncu- 
lus jb avidson & Humphre v9ll997l: Iz ethson 20011; 



. .. _ . . ^, Ishibashi et alj 

llOOaTlSmith et al..,2004.: .Nielsen et al.ll2007 a : .Mehner et al.ll2010h . 



Figure 12. Illustration of r] Car's binaiy orbit (inset, yellow) on the sky 
relative to the Homunculus nebula for a binary orientation with i = 138°, 
6 = 7°, and PAz = 317°, which lies near the center of our best-fit range 
of orbital parameters. The +z orbital axis (blue) is closely aligned with the 
Homunculus polar axis in 3-D. rje orbits clockwise on the sky relative to rfj^ 
(black arrows in inset), and apastron is on the observer's side of the system. 
The semi-major axis (+a;-axis, red) runs from NW to SE on the sky, while 
the semi-minor axis (-|-iy-axis, green) runs from SW to NE. North is up. 



way to absorb or diminish the extended, optically-thin red-shifted 
emission while allowing the blue-shifted emission component to 
reach the observer. The problem is that extended, red-shifted emis- 
sion is observed at other slit PAs, even for similar orbital phases. 
Therefore, the absorption/diminution would have to be restricted to 
very specific, extended regions on the sky, which seems unrealistic. 
Given all of these issues, as well as those concerning the observed 
phase dependence of the emission, an orientation that places apas- 
tron on our side of the system seems unavoidable. 



5.1.1 Obsen'ational Support from the Weigelt Blobs 

While a detailed review of rj Car's spectral variability is beyond the 
scope of this paper, it now seems well established that observed, 
phase-dependent, narrow (< 50 km s~^) high-ioniz ation forbidden 
emission lines form i n Weigelt blobs B, C, and D ((Davidson et al] 
1995l;IZethsonll200ll ; IZethsonet al]|201ll ; iDamineli et al.ll2008allbl: 
Mehner et al.ll2010h . [llumination of the Weigelt blobs by rjB is re- 
quired for the formation of these narrow lines, which are present 
during most of rj Car's 5.5-year orbit, and fade during perias- 
tron passage, grad ually returning to their 'normal' stren gth a few 
months afterward dVemer et al.ll2005hlMehner et al.ll201(]h . 

Our proposed orientation and direction for rj Car's binary or- 
bit are consistent with, and provide an explanation for, the time 
variab ility of the narrow-line emiss ion features seen in the Weigelt 
blobs. iDavidson et all ( [1991 1 19971) found that Weigelt blobs B, C, 
and D are located within 0.15" to 0.3" NW of the central stel- 
lar source near the Homunculus equatorial plane (see Figure |2ll. 



Since the Weigelt blobs are on the same side of the system as the 
observer, and because rye is required for the formation of the nar- 
row, high-ionization forbidden lines that form in the blobs during 
most of the orbit, one easily concludes that rjB must be on the same 
side of the system as the observer during most of the orbit. 

Moreover, the Weigelt blobs are located in the NW quadrant of 
the system. For our proposed orbital orientation, during most of the 
orbit, the low-density wind cavity created by rjB is open toward and 
pointing NW, in the direction of the blobs. The low-density cavity 
thus provides a path for the ionizing photons from rjB to reach the 
Weigelt blobs and produce the high-ionization emission. Over the 
course of the 5.5-year orbit, rjB gradually moves from SE to NW on 
the sky when going from periastron to apastron, and then from the 
NW back to the SE when moving from apastron to periastron. The 
orbital motion of tjb and gradual embedding in the wind of tja close 
to periastron leads to changes in the spatial extent and direction 
of the photoionization region, which causes observed variations in 
the narrow, high-excitation emission in the blobs. During periastron 
passage, the ionizing flux from tjb shuts off, and the high-excitation 
emission from the blobs fades. It takes several months after perias- 
tron for 77b to restore the large photoionization volume that faces 
NW, and thus the high-excitation emission in the blobs. 

While omitting some of the details, this qualitative picture is 
consistent with the known behavior of the Weigelt blobs. The pro- 
posed clo ckwise direction o f rjB's orbit is further supported by the 
images of ISmith et all | |2004) . which show excess UV emission to 
the SW of rj Car just before periastron, and to the NE just after. A 
complete review of all of the available observations of rj Car and 
how they support the derived orientation and direction of the binary 
orbit is obviously not possible here. However, we note that the ori- 
entation shown in Figure[T2]and clockwise motion on the sky of rjB 
are consistent with all known observations of rj Car to date. 



5.2 Implications for Theories for the Formation of the 
Homunculus Nebula and the Nature of t7a's Wind 

A variety of models, ranging from the interaction of stellar winds 
with differing speeds and densities, to binary interactions, and 
even a binary merger, have been proposed for explaining b oth the 
Great Eruption and the bipolar shape of the resulting nebula jSmithl 
I2OO9I) . Observational studies of the Homunculus have firmly estab- 
lished that its polar axis is orientated on the sky at a PA 312°, 
and that its inclination is 42° ± 1° l lDavidson e7ai]|200ll : TSmitbl 
I2006bl . l2009il . As shown in Figure [121 the orbital axis of the rj Car 
binary for our best-fit range of orientation parameters is closely 
aligned in 3-D with the Homunculus polar axis. Such an alignment 
has important implications for theories for the formation of the Ho- 
munculus and/or the present-day shape of 77a's wind. 

HST/STIS long-slit spectral observations of the Homunculus 
by Smith et al. (2003a) indicate that the reflected stellar spectrum 
over the poles has stronger and broader absorption in Ha, implying 
a denser, polar outflow. Very Large Telescope (VLT)A'ery Large 
Telescope Interferomet er (VLTI) observa t ions o f the optically thick 
stellar wind of rjA by Ivan Boekel et al.l | |2003|) and IWeigelt et al.l 
(2007) also indicate that ?7a's wind could be prolate in shape and 
aligned at the same PA on the sky as the Homunculus. Many have 
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interpreted these observational studies as evidence for tiie current- 
day wind of r/A being prolate, with a bipolar form and orientation 
similar to the Homunculus. Thus, one explanation offered for both 
the Homunculus' formation and the nature of the present-day wind 
is that r]A is a rapid rotator (Owocki & Gayley 1997; Maeder & 
Meynet 2000; Dwarkadas & Owocki 2002; Smith et al. 2003a; 
Owocki, Gayley, & Shaviv 2004); the effects of gravity darken- 
ing on radiation-driven wind outflows from a rapidly rotating star 
suggested as a way of explaining both the high polar densities and 
velocities inferred in jja's extended wind. 

Ho wev er, usi n g 2-D radiative transfer models, iGroh et al] 
l l2010ij) and lOrohl ( 1201 ll) find that the density structure of t^a's 
wind can be sufficiently disturbed by 77B, thus strongly affecting 
the observed UV spectrum, optical hydrogen lines, and mimick- 
ing the eff e cts of fast rotation in the interferometric observables. 
iGrohetaljfcoiOiJ) further show that even if t^a is a fast rotator, 
models of the interferometric data are not unique, and both prolate- 
and oblate-wind models can reproduce the interferometric observa- 
tions. These prolate- and oblate-wind models additionally suggest 
that the rotation axis of t^a would not be aligned with the Homuncu- 
lus polar axis. Further complicating the situation are the recent re- 
sults of[Mehner et al. (201 1), who find no evidence for higher wind 
velocities at high stellar latitudes in H5 P Cygni profiles obtained 
during r) Car's normal state (see their Section 5). Together, these 
new results challenge the idea that the present-day wind of r]A is 
latitudinally-dependent and prolate in shape. 

In this context, it is noteworthy that the results in Section 4 
match the observations so well, given that the winds from both stars 
are assumed to be spherical in the 3-D simulations. New 3-D simu- 
lations with a latitudinal-dependent wind for r/A, and comparison of 
the resulting spectro-images to the observations, may provide addi- 
tional clues as to the nature of t^a's present-day wind. It is impor- 
tant to keep in mind though that these new results do not necessar- 
ily mean that t^a's wind was not latitudinally-dependent in the past. 
Latitudinal-dependent mass-loss from a rapidly-rotating r/A could 
have played a role in the formation of the Homunculus. However, 
the argument for this possibility should not be made based on the 
above-mentioned observations of t^a's current-day wind. 

If t/a's rotation axis is aligned with the orbital axis, and the 
star is a rapid rotator, then the extended wind should be prolate 
and aligned at nearly the same PA as the Homunculus. However, 
the orbital axis and rotation axis of 77A do not have to be aligned. 
This is another common assumption, with important implications 
for theories proposed to explain r] Car's spectroscopi c events 
dPavidson & Humphrevsl 1 19971 ; ISmith et al.l |2003a; M ehner et al.l 
I2OIII) . Future observations using diagnostics that are not signifi- 
cantly affected by rys are needed to determine whether or not rjA 
is a rapid rotator, and whether its rotation axis is aligned with the 
orbital axis. 

Alignment of the orbital and Homunculus polar axes strongly 
suggests that binarity played a role in the Great Eruption, and pos- 
sibly the smaller eruption that later formed the Little Homunculus. 
A binary merger seems unlikely giv en the multiple known large 
eruptions of the system ( ISmithll2009l) . The most likely situation in- 
volves some sort of interaction between rjA and 77B at periastron. 
ISmith et al.l J2003ah suggested a scenario in which a rotating rjA 
loses mass due to increased angular momentum caused by an in- 
teraction with rjB- Models for the periastron-passage triggering of 
r) Car's massive eruptions, and the requirement of a binary to ex- 
plain the bipolar shape of the Homunculu s, have also been propose d 
by Soker and colleagu es ( Kashi & S oke r 20 1 d ; fs okeill2004 r2007h . 

More recently, ISmith & Frew! j2011 ) presented a revised his- 



torical light curve for 77 Car, showing that two 100-day peaks ob- 
served in 1838 and 1843 just before the Great Eruption coincide 
to within weeks of periastron, provided the orbital period then was 
shorter than the current period by ~ 5%. The beginning of rj Car's 
lesse r outburst in 1890 als o seems to have occurred around perias- 
tron jSmith & Fre"wll201 ih . Based on these findings and other con- 
siderations, ISmithl ( l201lh proposes that a stellar collision occurred 
at periastron before and during 77 Car's Great Eruption, with rjB 
plunging deeply into the bloated photosphere of rjA- Is close align- 
ment of the orbital axis and polar axis of the Homunculus evidence 
for such a scenario? Only detailed theoretical modeling can answer 
this question, but given the now apparent alignment of the two axes, 
binary interaction scenarios should be seriously considered as pos- 
sible explanations for rj Car's multiple, massive eruptions. 



5.3 Probing Changes in r] Car's Stellar and Wind 

Parameters via Long-Term Monitoring of tlie Forbidden 
Line Emission 

The rj Car system not only has a centuries-long history of vari- 
ation, including two major eruptio ns, it exhibits shorter-term 
variability with a 5.5-y e ar period ( ^ Davidson & Humphrevsl 1 19971 ; 
iHumphrevs et alj I2OO8I ; ISmith & Frew! I2OIII) . This implies that 
rj Car's numerous spectroscopic events are also related to periastron 
passage (Damineli et al. 2008a b). Contrary to many expectations, 
rj Car's 2009.0 spectroscopi c event was conside rably different than 
that of earlier cycles (Clh lMehner et al]|201lh . The observed X- 
ray minimum in 2009 was substantially shorter than the minima 
in 1998 and 2003, by approximately one month (CIO; Cll). The 
emission strength of Ha and various X-ray lines also appears to 
have decreased by factors of order tw o (Cll). Recent HST WFPC2 
observations bv lMehner et aljf201 ih show that the minimum in the 
UV was much deeper (again by about a factor of two) for the 2009 
event than for the 2003 event. While still the subject of debate, 
the behavior of [He II] A4686 might have been different in 2009 as 
well, with a second ' outburst' occurring ~ 30 days after the first 
dTeodoro et al.ll201 iLlMehner et al.ll201 ih . 

The cause of this sudden change in the behavior of rj Car's 
spectroscopic event is still poorly understood. One proposed expla- 
nation is that the mass-loss rate of rjA has suddenly de creased by 
at least a factor of two (CIO; Cll; iMehner et al.ll2oTTh . However, 
this is far from confirmed and the exact reasons for such a drop 
are unknown. The results in this paper suggest a possible way to 
observationally test the idea that 77a's mass-loss rate has recently 
changed. All of the HST/STIS observations modeled were taken in 
or before 2004. The model in this paper assumes a mass-loss rate of 
10"'^ M0 yr~^ for t^a, which is based on earlier observational and 
modeling studies (Davidson & Humphreys 1997; HOI; H06). The 
good match between the synthetic spectro-images and the observa- 
tions suggests that the mass-loss rate of 77A was ^ 10""^ M0 yr~^ 
at the time the observations were taken. 

The results in Section 4 though are strongly dependent on the 
mass-loss rate assumed for rjA and the ionizing flux of photons 
assumed for rjB- If the ionizing flux of photons from rjB remains 
constant, but the mass-loss rate of 77A drops by a factor of two or 
more, the size of the photoionization region should increase consid- 
erably. Therefore, the spatial extent and flux of the observed high- 
ionization forbidden lines should also drastically change. The phase 
dependence of the forbidden emission would likely differ, with the 
high-ionization emission vanishing at later phases (compared to 
earlier orbital cycles) when going into periastron, and reappearing 
at earlier phases afterward. Observed variations with slit PA should 
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be different too, possibly showing longer spatially-extended struc- 
tures, and even new components at different vios- Changes in the 
wind terminal velocity of rjA could be similarly investigated using 
the maximum observed vios of the spatially-extended emission. 

Multi-epoch observations coupled with improved 3-D radia- 
tive transfer modeling of the high-ionization forbidden line emis- 
sion would also help in determining if there is a significant change 
in TjB 's ionizing flux of photons, mass-loss rate, or wind terminal 
speed. Constraints on the ionizing flux of ph otons from 7?b could 
be comp ared to stellar mo dels for a range of O ^Martins et al.l2005l) 
and WR ICrowthej2007l) stars, allowing one to obtain a luminosity 
and temperature for tjb- 



6 SUMMARY AND CONCLUSIONS 

A major goal of this paper has been to use the 3-D dynamical 
model and available HST/STIS observations of high-ionization for- 
bidden lines to constrain the absolute 3-D orientation and direction 
of rj Car's binary orbit. The spatially-resolved spectroscopic obser- 
vations obtained with the HST/STIS provide the crucial spatial and 
velocity information needed to help accomplish this. Since the for- 
bidden lines do not have an absorption component, they provide 
a clear view throughout rj Car's extended interacting winds. Each 
high-ionization forbidden line arises only in regions where (1) pho- 
toionization by rjs can produce the required ion, (2) the density is 
near the line's critical density, and (3) the material is at the appro- 
priate temperature. Our results are therefore not as dependent on 
the many complicated details and effects encountered when mod- 
eling other forms of emission and/or absorption, such as X-rays or 
spectral features of Hi, He I, or other species (HOI; H06; Nielsen 
etal. 2007b; PU). 

Synthetic spectro-images of [Felll] emission line structures 
generated using 3-D SPH simulations of rj Car's binary WWC and 
radiative transfer codes were compared to the available HST/STIS 
observations for a variety of orbital phases and STIS slit PAs. The 
model spectro-images provide important details about the physical 
mechanisms responsible for the observed high-ionization forbidden 
line emission, as well as the location and orientation of the observed 
emitting structures. Below, we summarize our key conclusions and 
outline the direction of future work. 

(i) Large-scale (~ 1600 AU) 3-D SPH simulations of rj Car's 
binary colliding winds show that during periastron passage, the 
dense wind of rjA flows unimpeded in the direction of system apas- 
tron (Figure[6). Following periastron, rjB 's high- velocity wind col- 
lides with this primary wind material, creating dense shells that 
expand outward during rj Car's 5.5-year orbital cycle. These shells 
eventually fragment, forming a pair of dense, spatially-extended 
'arcs' of primary wind on the apastron side of the system. After 
several orbital cycles, the arcs drop in density and mix with the 
surrounding low-density wind material from rjB ■ 

(ii) During most of rj Car's orbit, far-UV radiation from rjB pho- 
toionizes a significant, spatially-extended region on the apastron 
side of the system, including portions of the dense arcs of primary 
wind formed following periastron (Figure |6]l. Portions of the cur- 
rent WWC region and primary wind just beyond it are also pho- 
toionized by rjB- The density and temperature of the material in 
these regions are ideal for producing [Felll] emission (Figure[6}. 

(iii) During periastron passage, rjB becomes enveloped in the 
dense wind of rjA- This significantly reduces the size of the pho- 
toionization region, to approximately a few semi-major axes in di- 
ameter Since the large photoionization region has 'collapsed', the 



spatially-extended (> 0.1") high-ionization forbidden line emis- 
sion vanishes during periastron passage (0 ~ 0.99 — 1.07). It does 
not reappear until tjb completes periastron passage and restores the 
extended photoionization region, which can take several months. 

(iv) Synthetic spectro-images of [Fe III] A4659 show that most 
of the [Felll] emissi on originates in the central 0.15", in agree- 
ment with fis ure 8 of iMehner et"ai] ilOldi . Our 3-D model shows 
that this emission forms in photoionized material near the cur- 
rent WWC zone, which is not spatially-resolved by HST/STIS. 
Spatially-extended (> 0.15") [Felll] emission arises in photoion- 
ized portions of the expanding arcs of primary wind formed just 
after periastron passage. 

(v) Only spectro-images generated for lines-of-sight angled — 
— 15° to +30° prograde of the semi-major axis on the apastron 
side of the system are able to match the observations taken at (j) — 
0.976, slit PA — 38°. However, there is an ambiguity in the orbital 
inchnation i, with models for i = 35° to 50°, PA^ = 272° to 332° 
and i = 130° to 145°, PA^ = 282° to 342° both able to produce 
entirely blue-shifted arcs (Figure|7]l. 

(vi) Model spectro-images of multi-phase observations obtained 
at slit PA — —28° break the above degeneracy in i, fully con- 
straining rj Car's orbital orientation parameters. Given the uncer- 
tainties in some of the stellar/wind parameters used in the 3-D mod- 
eling, our suggested best-fit range of orientation parameters for the 
rj Car binary system are i « 130° to 145°, 6 ~ —15° to +30°, 
PAz ~ 302° to 327°. Therefore, the orbital axis of the rj Car bi- 
nary system is closely aligned in 3-D with the Homunculus polar 
axis, with apastron on tiie observer's side of the system and rjB 
orbiting clockwise on the slcy relative to rjA (Figure [T2t . 

(vii) All other orbital orientations, including that in which rjB is 
in front of rjA at periastron (9 = 180°), are explicitly ruled out as 
they produce model spectro-images that lack entirely blue-shifted, 
spatially-extended arcs, containing instead significant amounts of 
unobserved, spatially-extended red-shifted emission at = 0.976, 
sfit PA = 38° (FigurelDl. 

(viii) The 3-D model predicts that the orientation on the sky of 
the blue-shifted component of [Fe III] A4659 should not vary much 
during rj Car's orbit, stretching from NE to SW. However, the spa- 
tial extent of the emission should grow larger as the system moves 
from periastron to apastron, and decrease in size close to perias- 
tron. The orientation and spatial extent of the red-shifted compo- 
nent should behave in a way similar to that of the blue-shifted com- 
ponent, but pointing NW on the sky. 

(ix) The apparent alignment of the orbital axis and Homunculus 
polar axis implies a link between binarity and rj Car's numerous 
massive eruptions. Binary interaction scenarios should be seriously 
considered as possible explanations for the Great Eruption and for- 
mation of the bipolar Homunculus, and possibly the smaller erup- 
tion in 1890 that formed the Little Homunculus. 

(x) Future detailed 3-D simulations and radiative transfer calcu- 
lations, together with high-resolution spatial mapping of the high- 
ionization forbidden line emission with HST/ST\S, have the poten- 
tial to further constrain the temperature and luminosity of rjB ■ The 
strong dependence of the forbidden line emission on the mass-loss 
rate of rjA and ionizing flux of photons from rjB may provide clues 
as to the nature of rj Car's long-term variability, specifically, the 
cycle-to-cycle variations of the spectroscopic events. 

The 3-D dynamical modeling in this paper and the available 
HST/STIS spectral observations have increased our understanding 
of the rj Car system. Much work remains, however. Efforts are un- 
derway to improve both the 3-D hydrodynamical modeling and 3-D 
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radiative transfer approach. New 3-D simulations that include ra- 
diative cooling, radiation-driven stellar winds, gravity, and effects 
like radiative inhibition and braking are on the horizon. Prelim- 
inary result s of detailed 3-D ra diative transfer calculations using 
Simplex ( |Paardekoope^l2010l) are promising as well. With these 
new tools at our disposal, it should be possible to further constrain 
the stellar, wind, and orbital parameters of r) Car, setting the stage 
for orbital modeling to determine the stellar masses. 
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